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ABSTRACT 

We explore the transmission spectrum of the Neptune-class exoplanet GJ 436b, including the possibility that 
its atmospheric opacity is dominated by a variety of nonequilibrium chemical products. We also validate our 
transmission code by demonstrating close agreement with analytic models that use only Rayleigh scattering or 
water vapor opacity. We find broad disagreement with radius variations predicted by another published model. 
For GJ 436b, the relative coolness of the planet's atmosphere, along with its implied high metallicity, may make 
it dissimilar in character compared to "hot Jupiters." Some recent observational and modeling efforts suggest low 
relative abundances of H2O and CH4 present in GJ 436b's atmosphere, compared to calculations from equilibrium 
chemistry. We include these characteristics in our models and examine the effects of absorption from methane- 
derived higher order hydrocarbons. To our knowledge, the effects of these nonequilibrium chemical products 
on the spectra of close-in giant planets has not previously been investigated. Significant absorption from HCN 
and C2H2 are found throughout the infrared, while C2H4 and C2H6 are less easily seen. We perform detailed 
simulations of JWST observations, including all likely noise sources, and find that we will be able to constrain 
chemical abundance regimes from this planet's transmission spectrum. For instance, the width of the features at 
1.5, 3.3, and 7 /im indicates the amount of HCN versus C2H2 present. The NIRSpec prism mode will be useful 
due to its large spectral range and the relatively large number of photo-electrons recorded per spectral resolution 
element. However, extremely bright host stars like GJ 436 may be better observed with a higher spectroscopic 
resolution mode in order to avoid detector saturation. We find that observations with the MIRI low resolution 
spectrograph should also have high signal-to-noise in the 5-10 /im range due to the brightness of the star and the 
relatively low spectral resolution (R ~ 100) of this mode. 

Subject headings: planetary systems; James Webb Space Telescope, radiative transfer; stars: GJ 436, HD 209458 



(N 
> 

in 

(N 
O 

o 



1. INTRODUCTION 

The rise of exoplanet characterization since the initial dis- 
covery of the transit of planet HD 209458 b a decade ago 
dCharbonneau et all l2000b iHenrv et all l2000h has been truly 
stunning. The initial steps in atmospheric characterization 
of hot Jupiters came after the realization that the transmis- 
sion spectra of transiting planets would be diagnostic of the 
temperature and chemica l mixing ratios in these atmosphere s 
(ISeaeer & Sasselovll200Ct lBrowrJl2001t iHubbard et all 1200 lh . 
As with perhaps all subfields of astronomy and planetary sci- 
ence, the return on the investment in these myriad observations 
is amplified when the data sets are compared to models that aim 
to simulate the conditions in the atmospheres of these planets. 

The model atmospheres, whether ID or 3D, aim to 
predict the temperature structure, chemical mixing ra- 
tios, and wavelength dependent opacity, as a function of 
height. A comparison of transmission spectra data with 
models can enable constraints o n the mixing ratios of 
atomic and molecula r absorbers (e . g.lCharbonneau et al] r2002; 
Fortnev et all 120031; ITinetti et all 120071: ISwain et al.l 12 008; 
Sing et all 120091; iM adhusudhan & Seagerl |2009t iDesert et al l 
2009; Bu rrows et alj|2010l) . Of course it then directly follows 
that constraints on these atmospheres are model dependent, as 
it is likely that the choices that one makes in constructing a 
model atmosphere affect the calculated spectrum. In this excit- 
ing field, the great difficulty in obtaining high signal-to-noise 
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observations, along with the somewhat unconstrained nature of 
current atmosphere models, can make interpre tation difficult. 
This w as a gain typified by the re cent works of IStevenson et al.l 
d2010b and Beaulieu et al] d2010l) . who, based on separate data 
reductions and model fits, disagree on the probable mixing ra- 
tio of methane in the atmosphere of GJ 436b, which we dicuss 
below. 

At the forefront of exoplanet characterization, detections of 
lower mass objects have expanded the reservoir of constraints 
to parameter space beyond hot Jupiters, broadening the scope of 
planetary characterization. The premier hot Neptune, GJ 436b 
a 22.6 M® planet orbiting an M2.5 star, was de tected via ra- 
dial veloci ty reflex motion bvlButler et all d2004l) . and later re- 
visited by iManess et al.l d2007l) . The first ph otometric detec- 
tion of transits were obtained by Gillon et al. ( 2007b), provid- 
ing the missing link needed to determ ine the mass and in turn 
the bulk density. iGillon et all d2007al) and lDeming et all d2007h 
refined sys tem parameters with a Sp itzer transit lightcurv e at 
8 urn, and iDeming et ail d2007l) and iDemorv et alTd2007l) re- 
ported a detection of the planet's secondary eclipse at 8 /zm as 
well. A recent attempt at transit characterization was lPont et al.l 
(2009), who probed the 1.4 fim water band using NICMOS on 
board the Hubble Space Telescope. This was the first attempt 
at a multi-wavelength transmission spectrum obtained for GJ 
436b. Ground base d efforts have yield e d H- and K-band ra - 
dius measurements dAlonso et al. 1 120081 ICaceres etalJ 12009). 
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in addition to radiu s measurements from EPOXI obtained by 
iBallard et all d2010h in the 0.35 - 1.0 u m range. Using IRAC 
on-board Spitzer, [Beaulieu e t al.l d2010l) very recently obtained 
transit depth measurements at 3.6, 4.5 and 8 /im. 

In o rder to study the dayside of the planet [Stevenson et al.l 
d2010l) obtained secondary eclipse measurements in six 
Spitzer bandpasses, from 3 .6 to 24 /im. They and 
Madhusudhan & Seager] d2010l) interpret their results as provid- 
ing evidence of thermochemical disequili brium in GJ 436b's 
daysi de atmo sphere. The fits by iMadhusu dhan & Seagerl 
(201 (J to the IStevenson et al] (120101) observations are quite in- 
teresting in the mixing ratios of CH4, H2O, CO, and CO2 that 
they require. They postulate a metal-rich atmosphere, in line 
with our understanding of Uranus and Neptune. However, com- 
pared to equilibrium chemistry, CH4 is strongly depleted, while 
H2O is also depleted, and CO and CO2 are strongly enhanced. 
Higher order hydrocarbon molecules, e.g., C2HX or HCN, were 
not considered in their fits, even though these molecules are 
the first products resultin g from methane destruction via pho - 
tolysis dMoses et alj|2005l) or reduction dZahnle et alj|2009al) . 
and are strong absorbers in the near and mid infrared. When 
suggesting that a low CH4 abundance may be due to the pho- 
tolysis and/or vertical mixing, these higher-ord er hydrocarbon 
molec ules should be included. Alte r nativel y, iBeaulieu et al.l 
d2010h interpret the IStevenson et alj d2010l) data (some re- 
reduced by their team), as well as primary transit data, as po- 
tentially indicating a methane-rich atmosphere with a temper- 
ature inversion. This only further strengthens the points that 
differences between data reduction methods and models clearly 
impacts conclusions. 

The coming of the James Webb Space Telescope will bring 
the next advancement in understanding exoplanet atmospheres 
and exoplanetary atmosphere modeling, as it will help to re- 
duce the under-constrained nature of current models. For a 
planet like GJ 436b where there are already signs from pho- 
tometry that the atmospheric chemistry is complex, the spec- 
tral capabilities of JWST will enhance our understanding, such 
as the identification of more complex chemical composition 
regimes and temperature structures. Understanding the integra- 
tion times and applicable wavelength bands needed to resolve 
spectral features of particular atmospheres is key for planning 
future science initiatives with JWST. 

GJ 436b is a particularly interesting planet because it is the 
transiting Neptune-class object that orbits the brightest parent 
star. This means it may long remain the best studied extrasolar 
Neptune-class object. Given the poor si gnal-to-noise of t rans- 
mission spectroscopy obtained with HST dPont et al.l2 009). and 
the uncertainties that arise from using a small number of wide 
photometric bands strung together as spectra, it may well fall to 
JWST to enable robust atmospheric characterization. 

Deriving accurate transit depths can be a complex task. Prob- 
lems with Spitzer IRAC primary transit observations include 
the issue of variable stellar fluxes between transit visits, due to 
starspots. In the most common situation where a planet does 
not occult spots directly, the unocculted parts of the stellar sur- 
face appear less bright, due to the spots. Compared to a stel- 
lar surface that lac ks spots, one would derive a deep er tran- 
sit depth (see, e.g. iPont etal] 120081: [Agol et al.ll2010l) . There 
is also the difficulty of extracting the transit depth with great 



accuracy, which yields different groups to find different tran- 
sit depths, with th e same data sets (e.g.. IBeaulieu et ail [2008: 
iDesert et al.l 2009). However, given that the interpretation of 
transmission spectroscopy is in principle less sensitive to the 
atmospheric pressure-temperature profile than day-side emis- 
sion spectroscopy, transmission spectroscopy may lead to the 
most robust determinations of atmospheric abundances. 

Here we examine the transmission spectrum signatures of 
a variety of atmospheric chemistries, including those from 
thermochemical equil ibrium calculations, those favored by 
IStevenson et alj d2010h . and those that include abundant higher 
order hydro carbons not considere d by these authors, guided by 
results from Zahn le et al.l d2009al) . We combine these models 
with a detailed simulation of the NIRSpec and MIRI Low Reso- 
lution Spectrograph instruments that will be aboard JWST, for a 
realistic simulation of what further knowledge we may gain for 
this important planet. In §2, we make simple comparisons be- 
tween our models, analytic relations, and other published work, 
in order to establish the robustness of derived atmospheric pa- 
rameters, and to validate our code. We find large differences 
between our work and that of G. Tinetti and collaborators. In 
§3, we describe model transmission spectra for GJ 436b while 
investigating various chemistries. In §4 we present JWST sim- 
ulations of GJ 436b transmission spectra, and discuss the possi- 
bilities for identifying important atmospheric signatures in the 
near future. In §5 we discuss implications of our transmission 
spectra models as well as address future endeavors for atmo- 
spheric characterization via theoretical work. 

2. TRANSMISSION MODEL: DESCRIPTION AND VALIDATION 

We model the transmission spectrum of planets using a de - 
scendent of the code first described in Hubbar det alj (|2001). 
Here we ignore the effects of refraction and a glow of photons 
aroun d the planet's l i mb du e to Rayleigh scattering, both of 
which iHubbard et alj d2001l) found to be negligi ble for close- 
in plan ets. Later works using this code included iFortnev et al.l 
d2003l) . which investigated simple two-dimensional models of 
the atmosphere of HD 209458b, but included one planet-wide 
P-T profile, w ith 2D changes in the atmospheric opacities. 
IFortnev! (2005) examined the possible effects of cloud opacity 
for the slant vie wing geometry ap propriate for transits. We re- 
fer the reader to lFortnev et al.l d2010l) . which described the code 
in some detail, and extended our treatment to 3D planetary at- 
mosphere models. 

In either ID or 3D, atmospheric P-T profiles are lain atop 
an opaque atmosphere at a reference pressure of either 1 or 10 
bar. The radius at this pressure level is adjusted to yield the 
best fit to observations. Along 1000 light ray paths through 
the atmosphere parallel to the star-planet-observer axis, the lo- 
cal atmospheric density and opacity are each typically sam- 
pled at 1000 points along each rajQ For absorption, the wave- 
length dependent cross-section is calculated based on contri- 
butions from a variety of atoms and molecules. The abun- 
dances can be based either on loc al chemical equil i brium 
at a given atmosp heric P-T point dLodders & Feglevl 120021 
2006; Lodders|2009j), or they can be arb itrary. The Rayleig h 
scattering cross-section is described in IFortnev et alj d2010l) . 
and in practice for our H2/He dominated atmospheres, we 
find a cross-section of 1.645e-24 cm 2 molecule" 1 at 450 nm, 
and scale by A -4 at other wavelengths. Here we define the 

6 Some phrases used in the original description of the code in Hubbard et al. 12001) may lead to the impression that the code replaces the true geometry of the 
atmosphere with a slab having the same column density. In the Hubbard et al. 12001) work this "slab" approximation was only done for the simulation of the weak 
glow of multiply Rayleigh-scattered photons. This approximation has never been made for ID or 3D calculations of the transmission spectrum. 
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where a and a are the wavelength dependent cross-section, 
and a reference cross-section, respectively, and A and X are 
the wavelength and a reference wavelength, respectively. Given 
these constraints, the planet's radius can be written as 



dR„ kT 

— — = a — = aH, 

dlnX pg 



(2) 



wavelength-dependent transit radius as the radius where the 
total slant optical depth reaches .56, following the results of 
ILecavelier des Etangs et al. (2008bfl who use a method simi- 
lar to our former one. 

2.1. Opacities 

Opacities generally used in the field of exoplane t atmo- 
sphe res are discussed in d et ail in [S harp & Burrows (2007) 
and IFreedman et all d2008l) . IFreedman etal J d2008h outline 
the opacities that we use in modeling the atmospheres of hot 
Jupiters, other Jupi ter-class planets, Nep tu ne-class planets, and 
brown dwarfs (e.g. Fortnev et al.| [2008a ,b; Sau mon et al.i r2006; 
ICushing et al.ll2008 ). We will not repeat the discussions in that 
paper, but we will touch on the issues of water opacity and 
higher-order hydrocarbons in turn. 

It is well estab l ished theoretically ( Burrows et al. 19971: 
Marl ev et al] 1 19991: ISeager & Sasselovl 120001: [Sudarsk v et al j 
2000l; lBarman et alj2001l) that water vapor opacity is the domi- 
nant infrared opacity source in warm giant planet atmospheres. 
There is also inescapable observational ev idence that this is 
true for brown dwarfs (e.g. Kirkpatrick 2005|), an d this clearly 
appears to be true for hot Jupiters as well (e.g. ISwain et alj 
2008; Grillmair et alj 12008). Comprehensive ab initio calcu- 
lations of line lists of hundreds of mil lions of lines for H2O 
have been tabulat ed by, for inst ance, Partridg e & Schwenkd 
d 1997b and bv lBarber et alj ( 120061) . Both tabulations are widely 
used. One of us (R. S. Freedman) has done extensive com- 
parisons of these two particular line lists at the tempera- 
tures of interest for planets and brown dwarf s, those below 
2500 K, and these differences are described in IFreedman et alj 
(2008]) at being "slight." In the detailed fits of M.S. Mar- 
ley and collabor ators to L- and T-type brown dwarfs, there is 
no hint that the iPartridge & Schwenkd (1 19971) database is in- 
sufficient to match the high signal-to-noise medium-resolution 
NIR and mid-I R spectra that have been achieved for sc ores 
of objects (e.g. Cushing et al. 2008} IStephens et al. 2009]). In 
FigureQ]we show calculated absorption cross-sections at 1500 
K and 1 mb a r. Cle arly the [Partridge & Schwenkd d 19971) and 
Barber et al. (2006) line lists are nearly identical in this pres- 
sure/temperature regime. If transmission spectra calculated by 
two different atmosphere codes differ (see 82.3I >, the choice be- 
tween these two water line lists cannot be an important con- 
tributing factor. 

2.2. Transmission Model Validation 

We can validate the predictions of the transmission spec- 
trum code by turning to previo us work. In particular, 
ILecavelier des Etangs et al. (2008a) have shown that the re- 
lation between absorption cross-section, mixing ratio, atmo- 
spheric temperature structure, and transit radius can be treated 
analytically. There is a particuarly straighforward relation 
for the wavelength-dependent transit radius for an atmosphere 
that obeys a few simple constraints. These constraints are an 
isothermal temperature structure, a constant gravitational accel- 
eration with height, and an opacity cross section a that varies 
as 

a = a (X/X o r, (1) 

7 In our published work he have used two different methods to calculate the "transit radius." In our earliest work, (e.g. Hubbard et al. 2001; Fortnev et al. 2003) 
we used the profile of the slant optical depth vs. radius, set atop an opaque circle, to generate synthetic images from which we calculated the amount of stellar flux 
blocked by the planet. We then defined a larger opaque circle whose cross-sectional area blocked this same amount of light as the planetary model — the radius of this 
larger circle yielded our transit ra dius. More recently we have chosen the tau=0.56 level as the transit radius, simply for computation ease, but we find fine agreement 
between the two methods, as did Lecavelier des Etangs et al. 1 2008b). Generally, for simple tests we find excellent agreement between our recent work and that of 
T. Barman Barman (e.g. 2007), who use a method similar to our current one, as well as with E. Miller-Ricci Kempton (e.g. Miller-Ricci et al. 2009). 



where R p is the transit radius, A is the wavelength, k is Boltz- 
mann's constant, T is the temperature, p is the mean molecular 
mass, g is the surface gravity, and H is the scale height. For a 
pure Rayleigh scattering atmosphere, a = -4. In Figure [2] we 
present an isothermal model at T = 1500K, g = 25 m s" 2 , and 
p = 2.32, with all opacity turned off, save Rayleigh scattering. 
We find a model slope that is within 1 % of the analytic relation 
from Eq. ©, which we regard as excellent agreement. 

2.3. Comparison With Other Work for Simple Models 

In 32. 11 w e described our implementat ion of the water vapor 
line list of IPartridge & Schwenkd dl997l) . As a further test of 
the transmission code we can isolate specific wavelength re- 
gions where water opacity closely obeys the a-relation from 
Eq. ([TJ. In particular, in Figure [3]2 we show the absorption 
cross-section vs. wavelength at 1500 K and 10 mbar. We have 
over-plotted fits for a in three spectral regions. The bluest and 
reddest wavelength ranges have a large negative slope, while 
the middle wavelength range has a positive slope. If our trans- 
mission spectrum model is working correctly, we should be able 
to match the transit radius slope of dR p /dln\ from Equation 
(O, for an atmosphere with a constant gravity and scale height, 
with water vapor being the only opacity source. We choose 
1500 K and the surface gravity of HD 209458b, 980 cm s" 2 . 
This model is plotted in Figure[3j>. One can readily see that over 
the three defined wavelength ranges that our model matches the 
analytic relation. 

Our choice of an isothermal HD 209458b-like model was 
based on models p resented by G . Tin etti and collaborators in 
a recent paper by iBeaulieu et alj (120091) . These authors used 
Spitzer IRAC observations to measure the transit depth in 4 
bandpasses from 3 to 10 pm. Their nice model fit, compared to 
the data (shown in their Figure 10), allowed the authors to as- 
sert that water vapor was the main absorber in th at atmosphere. 
This m ay well be true. However, as shown in iFortnev et alj 
d2010l) . our own HD 209458b models were not able to repro- 
d uce the large variatio n in absorption depths. As discussed 
in IFortnev et alj (j2010), we are generally unable to match the 
much larger variation in transit radius of the m odels of Tinetti 
and collaborators (e.g. lTinetti et a l. 2007, 2010). Although our 
two groups use different water opacity databases (Tinetii et 
al. use the BT2 list) it does not appear that can be a contributing 
factor. 

To help sort out t his issue, we b ecame interested in sim- 
ple tests. Figure 9 of IBeaulieu et alj (|2009) additionally shows 
transmission spectra for isothermal model atmospheres of HD 
209458b at 1500, 2000, and 2500 K, with opacity due only 
to water vapor (with a mixing ratio of 4.5 x 10~ 4 , the same 
va lue we use her e ). We compare our 1500 K model to that 
of IBeaulieu et alj d2009l) . as well as the analyti c relations, in 
Figure |4] (The model from lBeaulieu et all (120091) was obtained 
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using a data extraction software package.) The differences are 
large. Here we are able to match the analytic relations, while 
the model from lBeaulieu et afl d2009l) cannot. 

It is not immediately clear what causes these dramatic dif- 
ferences between the two models. The large differences remain 
at 2000 K and 2500 K as well. As shown in Figure Q] differ- 
ences in water opacity databases cannot be a culprit. Differ- 
ences in abundances can generally not be a reason either, as a 
higher (lower) water abundance would move the transit radius 
up (down) at all wavelengths. Also, that issue was eliminated 
for the simple test presented here. We are left in the position 
of identifying what we believe is a problem, with Tinetti et al. 
models, but we are not in a position to speculate as to its cause. 

We only dwell on this issue at length because it is the match 
of models to observations that allows for the identification of 
absorption features, and the determination of the mixing ratios 
of specific components. We are left to doubt the validity of 
the Tinetti transmission models, and the derivations of atmo- 
spheric abundances from some of the papers in which those 
models were used. We further wish to stress that at this time 
we have only compared to transmission spectrum models, and 
not to day-side emission spectrum models. Certainly the com- 
parison between models is an area in need of future work. An 
interesting avenue would be to model the t ransmission spec- 
trum of so lar system planets, su ch as Earth dPalle et alj [2009) 
or Saturn (Nicholson et al. 2006), which we will pursue in the 
near future. Our description of our methods complete, we can 
now turn our application to GJ 436b. 

3. APPLICATION TO GJ 436B 

GJ 436b is a re latively small planet with a bulk den sity 
similar to Neptune dGillon et al.ll2007bt iTorres et alj|2008l) . It 
may be mostly composed of fluid water, but a layer of H- 
He dominated atmosp here is clearly ne e ded to account for 
the observed radius (billon et alj l2007bt lAdams et al] [2008; 
Nettelma mTet alj |20 10). GJ 436b is one of the least-irradiated 
transiting planets, which makes its atmosphere cooler than 
many other well-studied planets. Based on thermochemi- 
cal equilibrium models, the low temperatures suggests the 
dominant c arbon-bearing molecule in the gaseous envelope 
is methane dSpiegel et aljhoiOt iMadhusudhan & Seageri feo 1 ; 



iLewis et al.ll201 



The recent IStevenson et al.l d2010l) secondary eclipse mea- 
surements from warm Spitzer, however, have been interpreted 
by these authors as suggesting otherwise. An atmosphere who's 
carbon chemistry is methane-dominate would tend to yield a 
small flux ratio in t he 3.6 nm band, with m ore flux in the 4.5 
^m band. However. IStevenson et alj d2~010) report a strong de- 
tection at 3.6 /im, and non-detection at 4.5 p,m. This could be 
indicative of extensive methane depletion, which would allow 
one to probe deeply, to hotter gas, in the 3.6 /im band, while a 
large mixing ratio for CO and CO2, which both absorb strongly 
in the 4.5 /im band, could suppress flux in this bandpass 
dStevenson et alj|2010t IMadhusudhan & Seageril2010l) . Large 
CO and CO2 ab undances have been shown to be indicators of 
high metallicitv dLodders & Feglevll20()2t IVisscher et a!l 2006; 
IZahnle et al.ll2009bl) . 

In the favored scenarios, the CO/CH4 and CO2/CH4 mixing 
ratios are enhanced due to vertical mixing from hotter, CH4- 
poor gas below, along with the photochemical destruction of 
CH4 by incident UV photons. If relatively abundant CH4 is 
indeed destroyed, this will give rise to a whole host of higher- 



order hydrocarbons, whic h is well understoo d for our solar sys- 
tem's giant planets (e.g. Mo ses et al. 1 120051) . Detailed chem- 
ical models, including photochemistry and vertical mixing , 
were applied to cool transiting planets by Zahnl e et alj d2009a). 
which predicted the formation of abundant C2H2, C2H4, C2H6, 
and HCN. More recently, the theory of a methane-poo r atmo - 
sphere for GJ 436b was challenged by iBeaulieu et alj (|2010), 
who, with a combination of modestly different eclipse depths 
at 3.6 and 4.5 /im, compute emission and transmission spectra 
models that allow abundant methane. 



3.1. Abundances and Opacities at the Terminator 

Since a goal of the work is to explore the prospects for broad- 
wavelength-coverage spectra with JWST, we choose to explore 
a diverse set of model atmospheres. Our choices are guided by 
predictions from equilibrium chemistry, nonequilibrium chem- 
istry, and fits to published spectra of the planet. Cross sections 
for the main molecules expected from thermochemical equilib- 
rium are shown in Figure [5] There are certainly a rich number 
of molecular bands, particularly at wavelengths blueward of 5 
/im, where NIRCam and NIRSpec will be sensitive. 

In Figure [6] we present GJ 436b model transmission spec- 
tra, using these cross sections. One model is for a 30 x solar 
metallicity atmosphere in thermochemical equilibrium (black), 
as well as two models (red and blue) with abundances taken 
directly from the "red model" and "blue model" in Table 2 of 
IStevenson et alj d2010h . These mixing ratios are somewhat de- 
pleted in water, strongly depleted in methane, and CO and CO2- 
rich compared to our equilibrium calculations. In all cases, we 
include the atomic sodium and potassium abundances derived 
from the equilibrium model. For the three chemistry cases, we 
use two different P—T profiles. In the upper panel, a hotter 
"dayside average" profile (named "2 pi") is used, which sim- 
ulates inefficient day night energy redistribution. The bottom 
panel uses a cooler plane t-wide average profile (named "4 pi") 
(e.g.. lFortnev et al1l2005l) . The two P—T profiles are plotted in 
Figure [7] These models (which we will refer to as models 'a' 
through 'f'), as well as models discussed in the following sec- 
tions are described in Table Q] At the current time we do not 
investigate time-variable temperature structure or abundances. 
The only published model o f the atmospheric dynamics of GJ 
436b show little variability dLewis et al.ll2010l) . and there is no 
observational evidence as yet. For the two P—T profile cases, 
we adjust the 10 bar radius of all models so that the radii align 
in the optical. The main chemical difference between the hot- 
ter and cooler profiles is that the hotter model yields a larger 
mixing ratio of Na and K, yielding stronger features in the opti- 
cal. A structural difference is that the larger scale height in the 
warmer model leads to modestly larger changes in radius as a 
function of wavelength. 

Compared to our eq uilibrium model, the two models from 
IStevenson etal] d201oi) (models 'c' through 'f in Table [B 
show considerably smaller radii in the near and mid-infrared, 
which is predominantly due to the smaller mixing ratio of wa- 
ter. Therefore, a clear probe of the water mixing ratio is the 
near-IR radius, compared to the optical radius, where water va- 
por is much less important. As expected the differences be- 
tween the equilibrium model and the red/blue models are great- 
est in wavelengths where methane and water are the dominant 
absorbers, and smallest where CO and CO2 are the dominant 
absorbers. 
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3.2. Nonequilibrium Chemical Products: Absorption Features 

IStevenson etall d2010t) and iMadhusudhan & Seagerl d2010l) 
suggested that there could be evidence for a very low CH4 mix- 
ing ratio in the atmosphere of GJ 436b, due to a combination of 
vertical mixing and photolysis of CH4. However, they did not 
investigate how the methane-derived nonequilibrium chemical 
produc ts may affect t he spec tra of the planet's atmosphere. Re- 
cently, IZahnle et al.l d2009al) have investigated nonequilibrium 
carbon for isothermal "warm Jupiter" atmospheres (800 < T < 
1200 K) at a range of metallicities. Th ese temperatures are sim- 
ilar to those suggested for GJ 436b. IZahnle et all d2009al) find 
that methane is sustained at higher regions in the atmosphere, 
water is more stable, and OH and H2 quickly combine to form 
H2O and H. This effectively increases the C to O ratio leading to 
increased abundances of molecules such as HCN, C2H2, C2H4, 
and C2H6. Given the uncertainties in modeling this chemistry, 
and the wide range of atmosphere models that are consistent 
with the GJ 436b data to date, our next aim is to explore trans- 
mission spectra with a range of nonequilibrium ch emical prod- 
ucts, guided by the results of Za hnle et al.l d2009al) . 

The cross sections for these first generation products of 
methane, which are strong absorbers in the infrared, can be 
found in Figure [8] Data for C2H2, C2H4 and C2H6, which are 
likely incomplete, are from the HIT RAN database, whi le that 
for HCN is from the calculations of Harri s et alj d2008l) . The 
cross-sections are similar in magnitude to those of the equi- 
librium chemistry products (Figure |5). HCN and C2H2 have 
prominent features at 1.5, 3.3, 7, and 13 /im, C2H4 has a signif- 
icant impact to absorption at 9.5 /im, but C2H6 has very little 
effect, with most of its opacity residing between 10 and 13 pxa.. 
In Figure|9]we explore transmission spectra including the opac- 
ities with these nonequilibrium products. We begin with model 
'e' in from Table Q] also shown as the red model from the top 
panel of Figure|6] This model uses a dayside average P—T pro- 
file with t he same best fit mixing r atios as the "red model" from 
Table 2 of lStevenson et alj d2010l) . 

In addition to the original parameters of model 'e', we in- 
clude absorption from HCN, C2H2, C2H4 and C2H6, with mix- 
ing ratios of 1 x 10" 4 , 1 x 10" 5 , 1 x 10" 3 and 1 x 10" 8 re- 
spectively, shown as the cyan model in Figure [9] (model 'g' 
from Table Q]). We de rive these chemical abundances based 
on IZahnle et al.l d2009al) . where they depict mixing ratios as a 
function of height and K zz (eddy diffusion coefficient) for an 
atmosphere at 1000 K. The mixing ratios chosen are similar to 
what is expected for an atmosphere that is vigorously mixing 
(K /y =W r ). 

It is worthwhile to explore some variations on this model. 
One variation removes the blanketing effects of hydrogen 
cyanide (HCN), shown in purple (model 'h' from Table []]). 
This model illustrates the difference in the width of the features 
at 3.3, 7 and 13 /im, where HCN and C2H2 have overlapping 
opacity. When the features are thick, they are dominated by 
HCN opacity. When we remove HCN, we see thinner features, 
indicating that C2H2 dominates over HCN. This characteris- 
tic may allow further constraints on the mixing ratios of cool 
transiting planet atmospheres. In green is a model with these 
nonequilibrium products absent at pressures below 10 mbar 
(model T from TableQ]). The green model illustrates the condi- 
tion of modest vertical mixing, with K 77 on the order of 1 6 cm 2 
s" 1 , similar to that favored by Madhusudhan & Seager (2010). 
In this case, nonequilibrium products are not stable higher in 
the atmosphere than the ~ 10 mbar pressure level dZahnle et al.l 



2009a). The main absorption features in the infrared will be 
weakened, as the green model confirms. The strength of these 
features could be constraints on vertical mixing and K zz . 

Absorption from these nonequilibrium products may mask 
features that would otherwise show an under-abundance of wa- 
ter in the transmission spectrum, as Figure [9] shows. The red 
model, which generally shows smaller radii in the near- and 
mid-IR than in the optical, instead shows significantly larger 
radii when the nonequilibrium products are introduced. Spec- 
troscopy, rather than photometry, will be key towards disen- 
tangling the effects of various molecules on the transmission 
spectrum. 

3.3. Comparison with Data 



in MgureliUj we explore model nts to recent bpitzer 1KAC 
data for 3.6, 4.5, and 8 /mi dBeaulieu et al.ll20iol). as well as 
data from EPOXI in the 0.35 - 1.0 /zm range (Ballar d et al 
120 Id) . HST NICMOS in the 1.1 - 1.9 Mm rang e dPontetal 



2009), ground based H-band dAlonso et alj|2008l) . and ground 
based K-band dCaceres e t al. 2009). We plot our models in high 
resolution, and as band-averages, where appropriate. Look- 
ing at the optical and near infrared, the best match is to the 
model that uses abundances from equilibrium chemistry, shown 
in grey (model 'a' from Table[TJ. 

In the mid-infrared, for the Spitzer IRAC data we see the 
same trend that we found in jFortnev et al] d2010t) : in com- 
parison to iBeaulieu et all d2008l) data for HP 189733b and to 
Bea ulieu et all d2009l) "data for HD 209458b, our models can- 
not match the larger amplitude of the features-the implied 
dramatic change in absorption depth as a function of wave- 
length. The models o f G. Tinetti and collaborators, used in the 
IBeaulieu etalld201 Oh paper, do fit the observation reasonably 
well. However, as discussed in §2.3, we find that the Tinetti et 
al. models overestimate the amplitude of absorption features. 

It is clear that our models do not ag ree with the large peak to 
trough variation of spectra required by Beaulieu et alj d2010h to 
fit their data . If the methods employed bv lBeaulieu et alj ([2008, 
120091 1201 01) are correct, and the error bars are not underesti- 
mated, the results for all of these planets imply that dramatic 
revisions to models of these atmospheres are needed. Given the 
uncertainties in the reduction of IRAC transit data, it may be up 
to JWST to confirm speculations about the molecules present at 
the terminator of GJ 436b. 

4. PROSPECTS FOR JWST FOR GJ 436 

We now evaluate the observability of the differences in the 
GJ 436b models by comparing them through the eyes of JWST. 
We have developed a code that simulates JWST spectra by com- 
puting the number of photons detected using a model of the 
host star, a transmission model of the planet, and estimates of 
the total efficiency (detected electrons per incident photon) at 
each wavelength for the various JWST dispersive spectroscopic 
modes. Noise is also modeled and added to the simulated spec- 
tra. 

The star GJ 436 is relatively bright over the 0.7-5 /im spec- 
tral region, and the transmission models predict that GJ 436b 
will have absorption features from many species over this wave- 
length range. Therefore we illustrate the model similarities and 
differences with simulations of JWST observations using the 
NIRSpec R = X/5X ^100 spectroscopic mode over this spec- 
tral range. The double-pass CaF2 prism used in this mode pro- 
vides spectroscopic resolution varying from R ~ 30 at A ~ 1.2 



6 



Shabram, et al. 



/mi to R > 200 at A > 4.3 ^m. We approximate this with a 
fourth order polynomial fit over the 0.7 - 5 /im range, and we 
assume that the prism has total transmission efficiency of 0.81 
after two passes. We estimate t hat the optical efficie ncy of NIR- 
Spec's 14 reflective surfaces dte Plate et a n i200l is approxi- 
mately 0.58 over A = 1 - 5 /mi, consistent with the values calcu- 
lated by the NIRSPec team (P. Jakobsen, private communica- 
tion 2003) and assumed by Deming et al. (2009) after remov- 
ing the grating blaze function. We adopt a quantum efficiency 
of 0.75 across the entire A = 1 - 5 /mi spectral range, consis- 
tent w ith the NIRSpec detector requirements (Rau scher et alj 
2007). The telescope is estimated to have total reflectivity of 
0.9 across this wavelength range. Total efficiency was modeled 
to decrease linearly by a factor of 2.0 as wavelength decreases 
from 1.0 to 0.7 /mi, driven mostly by reduction in reflectivity in 
the 14 reflective NIRSpec surfaces. 

Like lDeming et al] (120091) . we assume there will be no losses 
from the 1"6 wide entrance slit and that the only significant 
noise sources are photon noise and systematic noise due to 
small guiding errors during exposures. Photon noise is sim- 
ulated by adding Poisson noise appropriate for the number of 
detected photo-electrons in each resolution bin. We adopt the 
system atic noise value of 5 x 10~ 5 estimated by [Deming et al] 
(2009). We do assume that this noise is Gaussian in its distribu- 
tion although Deming et al. (2009) found that it was somewhat 
non-Gaussian. Even with high precision JWST instruments, 
we will suffer systematic noise at these modest but significant 
levels. In the "1 -(in-transit/star)" computation, any additional 
natural or instrumental noise occurring at frequencies greater 
than the inverse of the transit observation period will impact the 
extracted spectrum. This simulation program was coded in C, 
and it uses the public domain RANLIB package for simulating 
photon noise and Gaussian systematic noise. 

A high fidelity stellar model of the GJ 436 host star was not 
readily available, so we used a model of GJ 41 1 which has M2 
V spectral type, similar to the M2.5 V type of GJ 4 36. Using 
our simulation code, we re- binned the lKuruczl ( l2009l) R = 1000 
model of GJ 411 to the instrumental resolution of the JWST 
NIRSpec prism at each wavelength interval over A = 0.7-5 /jm. 
Next, our code computed the number of stellar photons from 
this binned flux, reducing it by the ratio of the squared model 
planet radius divided by the squared stellar radius (assumed to 
be 3.2 x 10 10 cm) at each wavelength. We used a distance of 
10.2 pc to GJ 436 and an integration time of 1800 s for these 
calculations. This integration time is ~ 33% s horter than the 
2740 s duration of the transit dPont et al]|2009l) . We used the 
resultant simulated in-transit spectrum and the simulated stellar 
spectrum of equal integration time to compute the absorption 
depth at each wavelength, 1 - (in-transit / star). This is plotted 
for GJ 436b models in Figure [TT] 

Models 'a', 'c', and 'e' (shown in black, blue, and red re- 
spectively) are plotted in the top panel of Figure [TT] These are 
JWST simulations of the models in the top panel of Figure [6] 
Models 'e', 'g', 'h', and T are plotted in the bottom panel in 
red, purple, cyan and green respectively, and are JWST simu- 
lations of the models in Figure [9] The absorption features are 
labeled here for clarity. In particular, in the top panel the dif- 
ferences between the water rich model (black) and water poor 
models (red and blue) are readily apparent, as is the strong CO2 
feature at 4.3 /im and CO feature at 4.5 /im. In the bottom 
panel, the absorption features due to nonequilibrium HCN and 
C2H-2 are clearly apparent, as is CH4 absorption from 3-4 /im. 



The prospects for detailed characterization of this planet, and 
others with JWST, is good. 

Given the high brightness of GJ 436b (K = 6.1 or K AB = 7.9 
mag), it is likely that its observation will require use of a detec- 
tor subarray that is smaller in the dispersion direction than the 
~ 35 pixel length of t he complete 0.6-5 /im R ~ 100 spec- 
trum (Tumlinson 2008). Therefore, we find that acquiring the 
entire spectrum shown at the signal-to-noise in Figure [TT] may 
require 2 or 3 transits. These observations may be best acquired 
in the higher resolution R = 1000 mode for stars as bright as GJ 
436. 

At A = 5 - 10 /im MIRI low resolution spectrograph (LRS; 
R ~ 100) observations were also simulated for these models in 
a similar fashion using details of LRS models and actual mea- 
sured performance. These simulations are shown in Figure [12] 
and are another independent way of looking at the effects of 
nonequilibrium chemistry in the atmosphere of GJ 436b. In 
the top panel (models 'a', 'c' and 'e' are shown as black, blue, 
and red respectively), water is the main opacity source in this 
wavelength range. The bottom panel (showing models 'e', 'g', 
'h', and T as red, purple, cyan, and green respectively) shows 
the clear distinction of HCN and C2H-2 between the different 
models shown at 7 /im. Absorption from C2H4 is shown at 9.5 
/tm as well. Observations using MIRI will be the only way to 
probe the 9.5 /tm feature created by the presence of C2H4. The 
simulations were also made for a total integration time of 30 
m in transit and 30 m on the star. The flux of GJ 436 is less 
than 1 Jy over this wavelength range, faint enough for its entire 
A = 5 - 10 /im spectrum to be acquired simultaneously, which is 
an advantage MIRI will have over NIRSpec observations that 
require multiple transits to obtain the full spectrum. 

The noise included in the spectral simulations and shown 
in Figures QT| and [12] are likely lower limits to the actual 
noise recorded in JWST spectra. Common-mode low frequency 
noises due to the observatory will likely be eliminated by the 
differential measurement of the star and planet. All data for the 
simulated spectra shown are obtained on time scales of 2 - 3 
eclipses, so any global variation in the star on longer time scales 
(i.e., subsequent transits or eclipses) will also be removed. Vari- 
ation of the planet between transits will be recorded, and this 
could, in principle, limit the usefulness of co-adding multi- 
ple spectra. Variations in the JWST instrumentation, the host 
star, or the planet on time scales shorter than 2-3 transit 
or eclipse events will appear in the data as systematic or ran- 
dom noise. Nevertheless, we expect that JWST should be 
able to obtain high signal-to-noise exoplanet spectra given that 
Spitzer was able to o btain signal-to-noise approaching 10 4 (e.g., 
Machale k et al]|20Toi and references therein) with its older de- 
tectors, about a factor of 2 less than we predict for the system- 
atic noise of JWST. 

5. DISCUSSION AND CONCLUSIONS 

We have explored possible transmission spectra of GJ 436b. 
We have investigated equilibrium chemistry cases and have 
compared them to mixing ratios derived to fit recent secondary 
eclipse measurements in the near infrared. The differences in 
these mixing ratio regimes have a significant effect when com- 
paring the corresponding transmission spectra. The lack of wa- 
ter in some cases will be easily detected by JWST assuming the 
spectrum is not dominated by absorption from nonequilibrium 
chemical products of methane. For this reason, we have made 
the first attempt to include opacity from higher order hydro- 
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carbons for close-in giant planet atmosphere models. Specif- 
ically, we have included absorption from HCN, C2H2, C2H4, 
and C2H6. The special condition that favors production of these 
molecules the most is the relatively low temperatures found in 
this planet. 

HST and Spitzer may be limited in the capabilit ies needed 
to ob serve this planet's transmission spectrum. The iPont et al.l 
(2009) near-IR spectrum from Hubble suffered from instru- 
mental systematics that could not be overcome. Reduction of 
Spitzer IRAC transit data suffers from the problem that dif- 
ferent groups reducing the same data sets achieve different re- 
sults, with errors b ars that do not overlap dBeaulieu et al.l l2008; 
Desert et al. 2009). Transit depths observed at different epochs 
allow for stellar vairability to complicate the interpretation. 
These observations are extremely challenging, such that many 
results for GJ 436b, and many of the transiting planets, must be 
regarded as provisional. 

We have shown that JWST should be able to make important 
breakthroughs in answering questions about the chemical mix- 
ing ratios in the planet's atmosphere. Specifically, the presence 
of higher order hydrocarbons in the transmission spectrum of 
this planet would reinforce the recent claim of non-detection of 
methane in the emission spectrum. We also may be able to con- 
strain the composition and structure at the planetary limb. If hy- 
drocarbon chemistry is present, it will be possible to constrain 
the abundances of these species. For example, HCN and C2H2 
have strong absorption overlapping at 1.5, 3.3, 7, and 13 pm. 
The thickness of the feature in each wavelength region will tell 
us which molecule is more abundant in the atmosphere. When 
there is minimal HCN, the signature of C2H2 in these wave- 
length regions will be characteristically thinner. This phenom- 
ena is depicted for JWST's NIRSpec in the lower panel of Fig- 
ure[TT]for the 1.5 and 3.3 /im features as well as for the MIRI 
LRS in the lower panel of Figure [12] at 7 pm. The strength 
of the features in the infrared produced by these hydrocarbons 
may also be instrumental in constraining information regarding 
atmospheric mixing. Weak features may indicate a small eddy 
diffusion coefficient (K zz ), pointing towards weak vertical mix- 
ing in the planet's atmosphere, with these molecules forming at 
higher pressure levels. 

In order to produce accurate model spectra of an exoplanet, 
including nonequilibrium chemical products (due to mixing and 
photochemistry) a self consistent code is needed. This would 
ensure that the molecular mixing ratios, and their impact on 
opacities, are consistent with the atmospheric P-T profile. This 
alone is difficult, but should be addressed in the future. In prin- 



ciple, all of this should be done in a three-dimensional model. 
Perhaps someday exquisite observations will warrant such a 
treatment. Clearly JWST will be an effective tool for reducing 
the under-constrained nature of current models, yet the bright- 
ness of GJ 436 may be a limitation. The NIRSpec prism mode 
will be able to observe a larger spectral range at high signal- 
to-noise with fewer transits for planets orbiting stars that are 
relatively dimmer than GJ 436. This issue however, does not 
pertain to obtaining spectra in the 5-10 pm range using the 
MIRI LRS. This study is the first in a series of papers that will 
investigate models of transmission and emission spectra of tran- 
siting planets, convolved with our realistic model of JWST ob- 
servations. For a complementary approac h over a large pha se 
space of JWST exoplanet observations, see Belu et al.1 (T2010). 

The model dependent nature of our understanding of the 
composition and structure of exoplanet atmospheres ensures 
that it will always be useful to refine atmosphere models. Thus 
developing, scrutinizing and refining these models brings us 
closer to the true nature of these planets. In §2 we validated our 
model against analytic relations, but found that a similar Tinetti 
et al. model could not match these same relations. Recent 
work on including nonequilibrium and pho tochemical prod- 
ucts dLiang et al.ll2004t IZahnle et all l2009bllah will yield bet- 
ter predictions of chemical mixing ratios. Transmission spec- 
trum models that incorporate the full 3D nature of the plane- 
tary atmos phere, which is particularly imporant at th e termina- 
tor region dFortnev et alj|2010t iBurrows et al.ll2010h will yield 
a more accurate understanding of the temperature structure and 
chemical abundances. The conditions of GJ 436b coupled with 
JWST will provide the chance to explore the possibility of a 
compositionally complicated atmosphere in great detail, ex- 
panding our knowledge of planets as unique to their particular 
circumstance. 
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Table 1 

Model Transmission Spectra for GJ 436b 
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1 x 10" 8 
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7 x 10" 


-4 


1 x 10" 7 


same as 


g, but absent at P < 10 mbar 
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FIG. 1 . — Absorption cross-section of water vapor at 1500 K and 1 mbar. The Partridge & Schwenke 1 1997) opacity is in red while Barber et al. 1 2006) in blue. 
Differences are very small. 
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FIG. 2. — Planet radius vs. wavelength for an isothermal, pure Rayleigh-scattering atmosphere. The analytical relation of Lecavelier des Etangs et al. (2008a) is 
shown as the solid line, while our constant-gravity isothermal model is the dashed line. The slopes of the lines agree to 1%. 
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FIG. 3. — (a) Cross-section of water vapor at 1500 K and 1 mbar, with analytic fits to a across three wavelength ranges. See Equations {TJ and {2) for details. 
Vertical dotted lines show the wavelength range of validity where we expect the analytic relation to hold, (b) The resulting planetary radius vs. ln(wavelength) of 
the model. Across the three wavelength ranges, the fit to the linear radius vs. ln( wavelength) is very good. Again, vertical dotted lines show the range of wavelength 
validity. 
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FIG. 4. — Absorption depth vs. wavelength for a 1500 K isothermal model, with the surface gravity of HD 209458b. In black is our model from Figure[3}>. The 
solid-colored curves are the analytic relations from Figure[5J> as well. In orange is a 1500 K, HD 20945 8b-gravity, model from Beaulieu et al. (2009). It is readily 
seen that our model presented here is a substantially better fit to the analytic relation. 
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FIG. 5. — Cross sections for the main equilibrium chemistry molecules at 10 mbar and 1500 K. NH3, CH4, FbO, CO and CO2 are represented by the black, red, 
blue, green and purple curves respectively. 
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FIG. 6. — Planet radius vs. wavelength for 2pi (dayside average) and 4pi (planet wide average) P—T profiles, at 3 different chemical abundance regimes. The 
black models are in chemical equilibrium with 30 X solar metallicity (models 'a' and 'b' from TablefTJ. The red and blue models adapt the abundances specified in 
Stevenson et al. (2010) (models 'c' through 'f'), and are best fits to secondary eclipse measurements at 3.6 and 4.5 /im. The Band average radii are plotted for these 
two bandpasses as squares in grey, cyan and orange corresponding to the black, blue and red models respectively. 
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FIG. 7. — Model pressure-temperature profiles for GJ 436b used in the current work. Profiles assume [M/H]=+1.5, or just above 30x solar. "2pi" is a dayside 
average profiles while "4pi" is a planetwide average. 



GJ 436b Transmission Spectra 



15 




0.6 



2 5 
Wavelength [/xm] 



10 



20 



FIG. 8. — Cross sections for the main higher-order hydrocarbon products at 10 mbar and 1500 K. HCN, C2H2, C2H4 and C2H6 are represented by the black, 
green, orange and cyan curves respectively. 
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FIG. 9. — Planet radius vs. wavelength using a 2pi (dayside average) P-T profile for the 'red' model described in Stevenson et al. (2010) shown in red (model 'e' 
from TablefT] Also shown are 'red' models with additional absorption due to nonequilibrium chemical products. In purple is model 'e' including absorption from 
HCN, C 2 H 2 , C 2 H 4 , and C 2 H 6 with abundances of 1 X 10~ 4 , 1 X 10~ 5 , 1 X 10~ 3 and 1 X 10~ 8 respectively (model *g' from TablefT). In cyan depicts this model 
again, removing absorption from HCN, revealing C2H2 as the dominant featu re at 3.3 fim (model 'h' from TablefT). In green we include chemical abundances that 
become absent above 10 mbar, a condition adapted from Zahnle et al. 1 2009a) (model T from TablefT). 
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FIG 10.— We present model fits to data f rom Ballard et al. 1 2010) (0.35 - 1.0 t an from EPOXI) IPont et al.l i2009T) (1.1 - 1.9 /im using HST MCMOSI. IAlonso et~aT] 
1 2008) (ground based H-band), Caceres et al] <2009l) (ground based K-band) and lBeaufieueFa l. 12010) (Spitzer IRAC 3.6, 4.5, and 8 fim bands). The grey model is 
a 30x solar metallicity model with a 2pi P—T profile (model 'a' from Table[T}. We show the red and cyan models from Figure[9]again here (models 'e' and 'g' from 
Table[T] respectively). It is clear that we are unable to reproduce the large peak to trough variation in spectra that allowed Beaulieu et al. 12010) to assert reasonable 
agreement of their models to current data. 




Wavelength [yu.m] 



FIG. 1 1 . — Simulated JWST NIRSpec prism mode observations of the absorption depth of the transiting planet GJ 436b relative to its host star at each wavelength. 
An integration time of 1800 s was used for both the in-transit and star-only spectra simulations. The models shown were previously plotted in Figures [6] and [9] 
Models 'a', 'c', and 'e' are shown in black, blue, and red respectively. Models 'e', 'g', 'h', and 'i' are shown in red, purple, cyan and green respectively. It is clear 
that we will be able to understand complex mixing ratio regimes with JWST. 
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FIG. 12. — Simulated JWST MIRI Low Resolution Spectrograph observations of the absorption depth of the transiting planet GJ 436b relative to its host star at 
each wavelength. Models 'a', 'c', and 'e' are shown in black, blue, and red respectively. Models 'e', 'g', 'h', and T are shownl in red, purple, cyan and green 
respectively.In the top panel, water is the main absorber in this wavelength range. In the bottom panel, HCN and C2H2 features are labeled at 7 fim, and C2H4 
absorption is labeled at 9.5 fim for clarity. 



